Adiabatic regularization of the graviton stress-energy tensor in de Sitter space-time 
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We study the renormalized energy-momentum tensor of gravitons in a de Sitter space-time. Af- 
ter canonically quantizing only the physical degrees of freedom, we adopt the standard adiabatic 
subtraction used for massless minimally coupled scalar fields as a regularization procedure and find 
that the energy density of gravitons in the E(3) invariant vacuum is proportional to H 4 , where H 
is the Hubble parameter, but with a positive sign. According to this result the scalar expansion 
rate, which is gauge invariant in de Sitter space-time, is increased by the fluctuations. This implies 
that gravitons may then add to conformally coupled matter in driving the Starobinsky model of 
inflation. 
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Interest in the back-reaction of particles produced by 
an external field has always been great in the last thirty 
years. This issue has been studied in cosmological space- 
times following the seminal works by L. Parker lj. The 
back-reaction problem for quantized test fields of dif- 
ferent spin has been examined in great detail over the 
decades for a de Sitter space-time 0- 

Gravitational back-reaction treated in a self-consistent 
way is much more difficult. Taking into account metric 
fluctuations together with matter fluctuations to first or- 
der is a textbook subject. To second order, the level of 
difficulty increases due to non-linearity and to the invari- 
ance under coordinate transformations. 

Vacuum space-times are the simplest arena in which to 
study the energy content carried by metric fluctuations, 
since the latter just reduce to gravitational waves (i.e. 
tensor modes) in the absence of dynamical matter. Even 
in such a simple setting, non-linear self-consistent solu- 
tions are quite surprising: the gravitational geon Q is 
one of the historic examples, which triggered more work 
in this area. 

Although gravitational waves are somewhat similar to 
massless minimally coupled test scalar fields, they are 
not exactly the same and the former may reveal quan- 
tum gravity aspects which are not proper of the latter: 
here we shall show this relevant difference in the de Sit- 
ter case. In this paper we study the energy-momentum 
tensor (EMT henceforth) of quantized gravitons and the 
scalar expansion rate in de Sitter space-time in a pertur- 
bative way. We shall consider only half of the de Sit- 
ter Carter-Penrose diagram, i. e. only the cosmological 
flat expanding branch. To our knowledge this case has 
been only approached by Ford 4], but not fully explored. 



Since inflation is very close to de Sitter, this calculation 
may also be useful in the context of back-reaction during 
inflation. For the case of inflation initiated by quantum 
anomalies of test fields Q, this would lead to the com- 
plete evaluation of the gravitational sector. Throughout 
units arc chosen such that h = c = 1 . 
According to the action 
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the Einstein equations are: 
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In the de Sitter space-time the only non- vanishing met- 
ric fluctuations are the dynamical degrees of freedom (the 
gravitons): 
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where a(t) = e Ht is the scale factor (H 2 = A/3), r\ is the 
conformal time and 7^ is the spatially flat metric (greek 
indices go from to 3, latin ones from 1 to 3 unless oth- 
erwise stated). The gravitons are traceless, transverse 
{h\ = , djh^ = 0) and therefore gauge- invariant with 
respect to tensorial spatial transformations. First order 
scalar metric fluctuations vanish in the absence of dy- 
namical matter, as already mentioned. 

In order to compute the graviton EMT we proceed as 
in textbooks 19: 
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ttG) 1 . The above expres- 
sion will become, after using the first order equations of 



where we have set M % — 



2 



motion: 
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where by the superscript (2) we mean terms which are 

(2) 

quadratic in the perturbation hij. R\tv can be found in 
Eq. (35.58b) of 0. 
We obtain 



From Eqs. JJjJ and l|10|) we see that the amplitudes h s ,k 
satisfy the same equation as massless minimally coupled 
scalar fields: 
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and the solution for the Fourier mode which becomes a 
plane wave for short wavelengths is: 
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This solution is valid for all values of k except for k = 
(j\0 (corresponding to a zero measure) in which case the 
solution is simply a space independent pure gauge. When 
averaged over the vacuum state annihilated by b the EMT 
of gravitons takes a perfect fluid form: 
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which is covariantly conserved in de Sitter space-time: 
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where ' denotes the derivative with respect to t. The 
action can also be expanded as S = S*' ' + S^- 2 \ where 
is the background value. The second order piece 
S( 2 \ omitting boundary terms, is: 
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Let us perform a Fourier expansion and consider only the 
physical degrees of freedom (polarization states) h + and 
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where e + and e x are the polarization tensors having the 
following properties (s = +, x): 
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These should be sufficient for our one-loop calculation on 
shell. Thus we do not concern ourselves with unphysical 
degrees of freedom and ghosts. This method of selecting 
only the physical degrees of freedom was used in first 
deriving the spectrum of gravitational waves and is 
used in computing the gravitational wave contribution 
to microwave anisotropies. 
On quantizing we have 
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where h s :k are the solutions given in Eq. (|15|) . The 
above expressions agree with the space averaged EMT of 
gravitons obtained in 8] . It is important to note that the 
term Hhh is reminiscent of a scalar field non-minimally 
coupled to gravity: however, a non-minimal scalar field 
would also have a mass term of order H (~ H 2 h 2 ) which 
instead is absent in Eqs. (|17I18|I . Owing to the presence 
of the term Hhh, the EMT of gravitational waves is not 
invariant under the transformation hij — > hij + const, 
which has generated so much activity in the context of 
massless minimally coupled scalar fields [13] ■ Of course, 
this is not surprising at all, since zero modes of gravitons 
are pure gauge, as already stated. 

In order to regularize bilinear quantities we proceed 
in the following way: we subtract the fourth order term 
of the adiabatic expansion for the solution to Eq. i|14|) 
with a mass term 0, . An EMT with a mass term 
(13) m which regularizes the adiabatic expansion is needed 
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in order to proceed with the adiabatic subtraction. It 
is then necessary to add a term to the energy density 
(e s ,k -* e s ,k(m) = e s , fe +m 2 |/i s fc| 2 /4) and to the pressure 
density (p s , k — > Ps ,k( m ) = Vs ,k + 5m 2 \h s fc 1 2 /12) in order 
to have a covariantly conserved EMT corresponding to 
()17I18JI with a mass term. The covariant conservation of 
such an EMT is given by 

e s M m ) + l e s A m ) + Ps A m )} = ( 19 ) 

and is equivalent to the equation of motion (|14J) with 
a mass term m 2 h s k ■ The regularized energy density is 
then given by 



which is only quadratically divergent in the ultraviolet 
(in contrast with the kinetic and gradient terms) and 
does not lead to any bad infrared behaviour (in contrast 
to |/i s fc| 2 ). The term in Eq. (|22[1 is therefore negative, 
with a sign which is opposite to that of the kinetic and 
gradient terms: this difference of sign is also reflected in 
the renormalized values, leading to the positivity of the 
energy density for gravitons. 

The result should be compared with that obtained for 
a massless minimally coupled scalar field in the Allen- 
Folacci (AF henceforth) vacuum 00: 
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where e^ k (jn) is the fourth order term of the adiabatic 
expansion 0. Let us note that, after performing the in- 
tegrals, we first take the limit m — > and finally remove 
the ultraviolet cut-off: in this manner the correct adia- 
batic subtraction is implemented. Indeed when this tech- 
nique is applied to a massless minimally cou pled scalar 
field all the results for an E(3) invariant-state [L| are re- 
produced. The above technique for computing integrals 
is therefore slightly different from the one we previously 
employed [T l ll5 | . 

For the EMT we finally obtain 



(£0) w hich corresponds to a contribution with a negative en- 
ergy density. On considering the AF vacuum for the cor- 
relator one obtains the following result 0, 0] : 
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This is our main result, which is in contrast with others 
claiming that gravitons decrease the effective cosmologi- 
cal constant at the one-loop order in de Sitter space-time 

SE3 

We stress that the result <|2"T|) is de Sitter invariant, al- 
though the vacuum chosen was E(3) invariant. The same 
happens for massless minimally coupled scalar fields ^3 • 
The interesting terms which break de Sitter invariance 
in the regularized value of EMT found in [Tol IT3| for 
an 0(4) state in the case of massless minimally coupled 
scalar fields are due to the use of closed spatial sections 
and are quickly redshifted after few e-folds of exponen- 
tial expansion [27|]. We also note that the result in Eq. 
(|21|l is obtained from the integration of the finite terms 
of the adiabatic expansion: this is also what happens for 
the averaged EMT of massless minimally coupled scalar 
fields. 

It is important to note that the renormalized EMT of 
gravitons in Eq. I|21|l has p = —p as equation of state. 
The relation p = —p/3 is the corresponding unrenormal- 
ized one for long- wavelength modes 

Let us again focus on the term 2Hhh: its Fourier com- 
ponent is 
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in which de Sitter invariance is broken in the standard 
way |23| • When the time derivative of Eq. (|24|l is con- 
sidered, one obtains 
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which is the same result as given by our method. Thus 
our approach is consistent with the choice of the AF vac- 
uum. The reason for the difference between l|21|) and 
(|23|l is the presence of the term 2Hhh for gravitons. Our 
mam result (EU can be easily verified by noting that the 
renormalized energy density of gravitons given by Eqs. 
(E|| is 
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The contribution of the term 2Hhh is positive and larger 
than the (negative) energy density of a massless mini- 
mally coupled scalar field. 

For the case of conformally invariant fields the EMT 
is independent of the vacuum state chosen and is fully 
given by the trace anomaly T: 
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with T given by |2lJ: 
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where C a pj5 is the Weyl tensor, which is zero for a met- 
ric which is conformal to Minkowsky as is Robertson- 
Walker (and therefore de Sitter). The coefficients a and f3 
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obtained by dimensional regularizations for scalar, four- 
component spinors, and gauge fields, are respectively 



a = < b , a = - < 11 . 
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For (massless) conformally coupled scalar fields one 
has: 

H 4 

(Tpu) ren = -g^v 9607r 2 ' ( 29 ) 

and therefore gravitons contribute, for example, as 361 
conformally coupled scalar fields! 

It is interesting to also investigate the effect of gravi- 
tons on the background space-time within the frame- 
work of second order perturbation theory for the Einstein 
equations, by evaluating the gauge invariant (for de Sitter 
space) geometric quantity O associated with the expan- 
sion rate of the universe (see for example |15|). Hence we 
consider the following second order metric fluctuations 
for the gauge fixed metric given in Eq. (|3*|): 

^ ) = y(^xf ) +^f + ^ ) ) ! 00) 

(2) 

where in the second line the vector \\ is divergenceless 

(2) 

and the tensor h) H is transverse and traceless. The ex- 
pansion scalar is defined by G = V m m m (where u M is a 
normalized vector field, u^u^ = — 1, defining the comov- 
ing frame) and simplifies to 

6 = ZH - ^h l} ti 3 - ZHa {2) + iv 2 /? (2) . (31) 

On using the Einstein equations one obtains the second 
order fluctuations as functions of the physical gravitons. 
In particular we find for the expansion scalar averaged 
over the vacuum: 

(1 pi tt2 \ 
1 + mo^5jj- < 32 > 

Thus we see that the choice of vacuum for the physical 
gravitons, which led to the result Eq. Q21|l correspond- 
ing to a positive cosmological constant contribution, also 
leads to a contribution of the same sign to the scalar 
expansion rate 0. 

Our result is only in apparent contradiction with the 
possibility that scalar fluctuations act against the acceler- 
ated expansion in chaotic inflation [l5j . A possible expla- 
nation of the difference is the stability of the space-time 



backgrounds: the inflationary era in scalar field driven 
universes is a transient state (local attractor), while de 
Sitter is a global solution. 

To conclude we have computed the regularized gravi- 
ton EMT in de Sitter space-time by quantizing only the 
physical degrees of freedom. We have found that the 
(one-loop) graviton contribution to the cosmological con- 
stant in the E(3) invariant-vacuum is positive, in contrast 
with that of a massless minimally coupled scalar field. 
This effect also appears in a second order perturbative 
analysis of the geometrical quantity O, which shows an 
increased expansion rate. According to this result, gravi- 
tons may then add to the trace anomaly of conformally 
coupled matter in driving the Starobinsky model of infla- 
tion The contribution of gravitons to the cosmologi- 
cal constant is not negligible and corresponds to a large 
number (361) of conformally coupled scalar fields. This 
contribution may also alter the inflationary phase of the 
Starobinsky model 0, |2^| since gravitons are not con- 
formally coupled, thus the back-reaction may alter the 
evolution of the gravitons themselves. 

Since gravitons are not conformally coupled [24j the 
averaged EMT may be state dependent. One may 
also worry about the problems concerning zero modes 
which p la g ued massless minimally coupled scalar fields 
[ToL ll'l flTl ITsf . Gravitational waves do not have zero 
modes (insofar as these correspond to a pure gauge), i. 
c. k > 0. Further all the contributions from the infrared 
to renormalizcd bilinear quantities which wc compute in 
this paper are finite (we may even include the contribu- 
tion from the k = mode since the measure of this point 
in Fourier space is zero). 

For a massless minimally coupled test scalar field in 
de Sitter space-time, it has been shown that the EMT 
evaluated in the AF vacuum is an asymptotic attractor 
among all possible vacua |25j | . If this were true for gravi- 
tational waves also, the result obtained for the E(3) state 
would be completely general and lead to an asymptotic 
value for the generalized anomaly j2f| Q 2 = 361/180 for 
gravitational waves. 

From the theoretical point of view it would also be in- 
teresting to see if our result (classical and quantum grav- 
ity in de Sitter space) can be related to conformal field 
theory, as suggested by the dS / CFT correspondence (2(| . 

Last, but not least, should the same result persists 
for cosmologies with H ^ and non-vacuum states for 
modes on large scales, it would be interesting to com- 
pute the contribution of cosmological perturbations to 
the present energy density. 
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